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Conventional spin-singlet Cooper pairs convert into spin-triplet pairs in ferromagnetic Josephson junctions 
in which the superconductor/ferromagnet interfaces (S/F) are magnetically inhomogeneous. Although 
much of the theoretical work describing this triplet proximity effect has considered ideal junctions with 
magnetic domain walls (DW) at the interfaces, in practice it is not easily possible to isolate a DW and 
propagate a supercurrent through it. The rare- earth magnet Gd can form a field- tuneable in-plane Bloch 
DW if grown between non-co -linearly aligned ferromagnets. Here we report supercurrents through 
magnetic Ni-Gd-Ni nanopillars: by field annealing at room temperature, we are able to modifj^ the low 
temperature DW- state in Gd and this result has a striking effect on the junction supercurrent at 4.2 K. We 
argue that this result can only be explained in terms of the interconversion of triplet and singlet pairs, the 
efficiency of which depends on the magnetic helicity of the structure. 



upercurrents in ferromagnetic Josephson junctions are strongly suppressed unless the superconductor/ 
ferromagnet (S/F) interfaces are magnetically inhomogeneous so that conventional singlet Cooper pairs can 
convert into spin-aligned triplet pairsand vice-versa^'^. Evidence for triplet- induced superconductivity in S- 
F hybrids was first reported more than ten years ago^ '*, predating the experimental proof of Josephson singlet Pi- 
coupling in S-F-S Josephson junctions (see, e.g., Ref. 5-17). A few years later, supercurrents in Josephson 
junctions with half-metallic ferromagnetic Cr02 barriers were also reported^^. These pioneering experiments 
had relatively low reproducibility because the S/F interfaces relied on the spontaneous appearance of inhomo- 
geneous magnetic states to generate triplet pairs. Recently, however, a series of experiments have been reported 
that demonstrate triplet pair supercurrent creation in S-F-S junctions by engineering the S/F interfaces to be 
magnetically inhomogeneous See also related works reporting evidence of induced triplet superconductivity 
using scanning tunnelling spectroscopy on S-F bilayers^^'^^, conductance measurements and Andreev spectro- 
scopy in S-F junctions^^"^\ and critical temperature measurements of S-F spin-valves^^. 

Due to the intrinsic magnetic inhomogeneity, a magnetic domain wall (DW) at an S/F interface should 
theoretically be an ideal triplet pair generator (see, e.g., Refs. 1 and 33-35). In most ferromagnets, however, 
the DW density even at the coercive field is not large enough to have a significant effect on the Josephson current 
of an S-F-S junction. This means that it is practically impossible to probe the effect that DWs have on the 
supercurrents that pass through them. As an alternative, intrinsically DW-like ferromagnets such as Ho can 
be used^^ and indeed our group^°'^^ (see also Ref. 38) and Sosnin et aU'^ have demonstrated results consistent with 
a triplet proximity effect in Ho-based junctions. Nevertheless, due to its high intrinsic anisotropy and interlay er 
exchange coupling energies it is not possible to unwind the helimagnetic structure without applying large 
(typically larger than 4 T) magnetic fields during transport measurements that would strongly suppress the 
superconductivity and any Josephson current in the superconductors studied; this makes it impossible to study 
directly how the magnetic structure controls the resultant supercurrent. A proper test for the role of DWs in 
mediating supercurrents requires therefore a magnetic structure in which the DW can be modified through 
magnetic field history. 

Exchange coupling at the interface between rare earth (RE) and transition metal (TM) ferromagnets results in 
an equilibrium antiparallel alignment of their moments. Non -parallel alignment of the TM layer moments in a 
TM-RE-TM trilayer can result in an in-plane Bloch DW in the RE^^ which can be modified by changing the 
alignment of the TM layers. Here we report supercurrents in Nb-Ni-Gd-Ni-Nb nanopillar Josephson junctions: 
depending on the Gd thickness, magnetic field annealing can significantly reduce the net barrier moment and 
substantially increase the junction supercurrent. 
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Figure 1 | Field-cooled temperature dependence ofthe magnetic moment 
per unit area for an unpatterned Nb(300 nm)-Ni(5 nm)-Gd(5.2 nm)- 
Ni(2 nm)-Nb(300 nm) film. The magnetic fields were applied in-plane as 
shown by the illustration. The inset shows the magnetic moment per unit 
area versus in-plane magnetic field for the same film at 9.8 K and 
comparative data from a Nb-Ni-Nb film at 100 K with a Ni thickness of 
7 nm. Solid lines are a guide to the eye. 

It has been shown previously that for barriers consisting of two 
ferromagnetic layers, parallel and antiparallel alignments of the layer 
moments result in different critical currents^^''*^ Since the magnetism 
in that case is always co-linear (i.e., there is no inhomogeneity of the 
form necessary to initiate spin-aligned triplet pairing), the larger 
critical current observed for antiparallel alignment originates from 
the "unwinding" of the singlet pair dephasing. We will show later 
that this effect cannot explain the observations reported here. 

Results 

Magnetic properties. For our experiment, a series of samples were 
prepared with different Gd thicknesses (dcd) between 1.1 nm and 
5.2 nm. The top and bottom Nb layers were both 300-nm-thick, 
and the top and bottom Ni layers had thicknesses 2 nm and 5 nm, 
respectively. Figure 1 shows the field-cooled temperature dependence 
of the magnetic moment per unit area for an unpatterned Nb-Ni-Gd- 
Ni-Nb film with a 5. 2- nm- thick Gd layer, measured using a vibrating 
sample magnetometer in which the field was applied in-plane as 
shown. The data shows a compensation point near 50 K where the 
moments of the antiparallel- aligned Ni and Gd are equal. Results for 
other Gd thicknesses are similar. A similar compensation point was 
reported by Barth et al for Ni-Gd bilayers^^. 

Nanopillar Josephson junction characteristics. Nanopillar Nb-Ni- 
Gd-Ni-Nb Josephson junctions were fabricated using a focused ion 
beam technique (see Methods Section and Refs. 43 and 44 for further 
information). Electrical measurements on junctions in a four-point 
current-biased configuration were performed using a magnetically- 
shielded dip-stick probe in a liquid He dewar. Current-voltage (IV) 
curves were measured with a lock-in amplifier from which the nor- 
mal state resistance (Rn) of a junction was measured at high voltage 
with a well-defined noise error. 

To confirm a Josephson effect in these junctions, we measured the 
effect of an in-plane magnetic field (H) on a junction's critical current 
(Ic). Fraunhofer-like oscillations of the critical current were obtained 
in all junctions, as shown in Fig. 2. The Ic(H) curves are hysteretic in 
magnetic field. This is because the total flux within the junction is the 
sum of that due to the applied field and the intrinsic flux generated by 



the net magnetisation (Mx) of the composite Ni-Gd-Ni barrier. For 
the simple case in which M^ follows a square hysteresis loop, the 
relative field-shift (AH) in Ic(H) is given by aM^, where a is a 
geometric flux- coupling factor. In other words, the maximum critical 
current (Ijviax) is offset by a finite field ± AH that is proportional to 
the net magnetisation of the barrier. 

The Ic(H) patterns of all junctions were measured in three stages. 
Stage 1 was the virgin state in which no magnetic field was applied to 
the junction before cooling down to 4.2 K [Fig. 2(i)]. For stage 2 the 
junctions were warmed to room temperature and a 0.1 T magnetic 
field was applied in-plane. After switching off the field, the junctions 
were cooled to 4.2 K and Ic(H) was re-measured [Fig. 2(ii)]. Finally 
(Stage 3), the junctions were warmed to room temperature and a 
1.0 T magnetic field was applied in-plane. After switching off the 
field, the junctions were cooled to 4.2 K and Ic(H) was re-measured 
[Fig. 2(iii)]. 

From Fig. 2 it can be seen that, with exception of the dcd =1.1 nm 
data, AH and iMax in the Fraunhofer patterns are dramatically chan- 
ged by the room temperature field treatment where, for progressively 
larger applied fields, iMax increases and AH decreases. This is more 
clearly shown in Fig. 3(a) - 3(d) where we plot the relation between 
Ijviax and AH for a number of different devices for each Gd thickness 
(to compare junctions with different cross-sectional areas and Gd 
thicknesses, Ijviax was multiplied by Rn to give the maximum char- 
acteristic voltage IjviaxRN as a function of AH). In Fig. 3(e) we show 
the dependence of the average maximum and minimum values of 
iMaxRN as a function of Gd thickness. 

In nanopillar devices in which the barrier is a single ferromagnetic 
layer, the value of AH implies that the barrier is single domain^°. For 
the Ni-Gd-Ni devices reported here, AH is dependent on magnetic 
history and so must have a more complicated magnetic structure 
reflecting a competition between the exchange coupling within the 
Ni and Gd layers, the interlayer exchange coupling, the magneto- 
static coupling between the layers, and the magnetic history of the 
samples. In the unpatterned films, the magnetostatic coupling is 
small compared to the various exchange energies and at low tem- 
peratures these exchange energies are minimised when the angle 
between the Ni layers (0) is zero and the Gd is homogenously mag- 
netised in the opposite direction. The compensation point in the data 
of Fig. 1 demonstrates this; henceforth this magnetic structure is 
referred to as the co -linear ferrimagnetic configuration. In the pat- 
terned nanopillars, the magnetostatic energy associated with the 
fringing fields is significant^^: at high temperatures where the Gd is 
non-magnetic, this favors 0 = 180°. However, at low temperatures, 
the cross-coupling between the Ni and Gd layers adds an additional 
exchange energy which acts to destabilise this arrangement. 

Here we focus on the formation of an in-plane Bloch DW in Gd as a 
means of reducing the magnetostatic energy [see Fig. 4]. It was shown 
earlier^^ that the energy of a DW in such structures is modified by the 
antiferromagnetic exchange interaction between transition-metals and 
Gd. Experimentally it was demonstrated that in micrometre- scale Py- 
Gd-Py (permalloy, Ni8oFe2o) trilayer mesas the magnetostatic energy 
is eliminated through the formation of a 180° in-plane DW within the 
Gd layer even for Gd thicknesses as small as 2 nm, which enabled the 
100-nm-thick Py layers to align antiparalleP^. The Ni layers are much 
thinner in our devices than in those in Ref 39, but the nanopillar 
dimensions are correspondingly smaller and so the magnetostatic 
energy remains substantial relative to the exchange energy. This sug- 
gests that, with increasing Gd thickness, there should be a crossover 
from the co -linear ferrimagnetic state to one in which a Gd DW is 
nucleated to reduce the magnetostatic energy. 

We now present a general expression for AH (dcd) in a Ni-Gd-Ni 
nanopillar barrier with and without an in-plane DW in Gd. To ac- 
commodate the misalignment between the Ni layers, we assume that 
the Gd magnetisation is perfectly helical, i.e. the in-plane magnetisa- 
tion angle (p(z) varies linearly along the junction axis z (orthogonal to 
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Figure 2 | Critical current modulation of Nb-Ni-Gd-Ni-Nb junctions, (a-c) Fraunhofer patterns at 4.2 K for samples with different Gd thicknesses 
(dcd) at stages 1,2, and 3, as described in the main text. The magnetic field was applied parallel to x [see inset in (b)(i)] except in (b)(ii) where the field wa; 
applied parallel to y. 
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Figure 3 | Supercurrent dependence of Nb-Ni-Gd-Ni-Nb junctions at 4.2 K. (a-d) The dependence of the Imex^n product on the net intrinsic field of 
the Ni-Gd-Ni barriers for samples with different Gd spacer thickness; (e) The average minimum and average maximum values of ImexRn versus Gd 
thickness measured at stages 1 and 3 respectively; the arrowed squares represent thicknesses at which no supercurrent was measurable. 
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Figure 4 | Possible low temperature magnetic configurations of Nb-Ni- 
Gd-Ni-Nb nanopillar junctions, (a) A co-linear ferrimagnetic Ni-Gd-Ni 
trilayer structure, (b) Antiparallel Ni layers with Gd forming an in-plane 
180° Bloch domain wall, (c) Parallel Ni layers with Gd forming an in-plane 
360° Bloch domain wall, (d) An in-plane flux-closed configuration. 



the layer interfaces). The shift, AH, in the Fraunhofer pattern is 
proportional to the net moment which is then given by: 
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where y = Mod / ^m- From the magnetisation data in Fig. 1, we 
estimate y to be approximately 10. 

In Fig. 5, we have plotted theoretical values of AH(dGd) on Gd 
thickness from Eq. (1) for three different values of 0, together with 
the average maximum and minimum values of AH obtained from the 
Ic(H) measurements. The constant of proportionality in Eq. (1) is 



0°Co-linear 
ferrimagnetic DW> 



80° DW 




Max. I AH I • 

Min.|AH| o 



^Gd ^^^^ 

Figure 5 | Intrinsic field of nanopilUar Nb-Ni-Gd-Ni-Nb junctions at 

4.2 K. The solid curves show the theoretical dependence of AH from Eq. 
( 1 ) on Gd thickness for different Bloch domain wall (DW) angles. The 
filled and hollow data represent the maximum and minimum values of AH 
obtained from the Ic(H) measurements shown in Fig. 2 while the shaded 
region represents the spread in AH. (Inset) The theoretical dependence of 
AH on DW angle for different Gd thicknesses. 



chosen in such a way that the co-linear ferrimagnetic curve coincides 
with the data point at dca = 0. No other fitting parameter is used. 

The CO -linear ferrimagnetic curve passes near the dcd — Tl nm 
data point, which is in agreement with our initial assumption that for 
small Gd thicknesses it is energetically unfavourable to nucleate a 
DW. For dcd > Tl nm, the largest measured values of AH are closer 
to the 180° DW curve than the co -linear ferrimagnetic curve, imply- 
ing the formation of a DW for these samples. Following high field 
annealing, the data points for dcd > Tl nm He close to the curve 
corresponding to a 450° DW. As shown in the inset to Fig. 5, the 
minimum barrier flux occurs around this angle and so this config- 
uration has a minimum magnetostatic energy. Clearly this structure 
has a much larger magnetic inhomogeneity so this, within a triplet 
picture, is consistent with the enhanced critical current. 

We know of no way of experimentally confirming these magnetic 
configurations in isolated nanopillar devices. However, we believe 
that sufficient supporting information exists to make the DW scen- 
ario plausible. Firstly, it has been confirmed both magnetically and 
electrically that a DW can be nucleated in Gd layer as thin as 2 nm^^. 
Secondly, the fact that the barrier magnetism can be changed by 
magnetic annealing implies that the Ni magnetisation is irreversible 
at room temperature (RT) even in the absence of a significant 
exchange interaction with the Gd: applying progressively larger fields 
must make this RT state closer to 0 = 0, while the virgin state is likely 
to be closer to 0 = 180° to minimise the magnetostatic energy. 
Finally, on cooling through the Curie temperature of Gd, various 
DW angles can be nucleated from a given starting value of 0: for 
example, for 0 = 0, either a co -linear ferrimagnetic configuration or a 
360° DW could form. As shown in the inset to Fig. 5, the latter has a 
much lower moment (and hence magnetostatic energy) for larger Gd 
thicknesses, and the moment can be further reduced by a relaxation 
to a DW angle greater than 360°. 

Discussion 

The analysis above provides a direct link between the rotation angle 
of the internal DW (as measured by AH) and the Ic- Before we can 
draw any conclusions regarding the role of unconventional super- 
conductivity in the behaviour, it is important to eliminate a rather 
trivial effect of flux in controlling the critical current via the standard 
Josephson relation Ij = Iq sin cp, where 0 is the local phase- difference 
across the barrier. For a square magnetic hysteresis loop, a net barrier 
flux gives rise to a standard Fraunhofer pattern with an offset AH 
which reverses at the coercive field. The important aspect of this is 
that the maximum critical current corresponding to the complete 
cancellation of the barrier flux by the magnetic field is identical to Iq. 

This cancellation is possible only because the barrier magnetisa- 
tion is laterally homogeneous. For inhomogeneous magnetisation, 
such as in -plane flux closure structures of the type illustrated in 
Fig. 4(d), the applied field cannot eliminate a phase variation across 
the junction and hence Imsx < loJ this has been experimentally 
demonstrated in much larger junctions than ours^^. This argument 
can be extended to barriers consisting of exchange- coupled magnetic 
multilayers of the type considered here. The conclusion is that the 
lateral magnetic inhomogeneity associated with the in-plane flux 
closure results in a suppression of the maximum critical current with 
respect to the non-magnetic or single-domain magnetic states: in 
other words, a dependence of iMax on AH which is opposite to that 
observed here. 

We now consider the configuration which we have proposed for 
our devices, in which an in-plane DW lies in the Gd layer [see 
Figs. 4(b) and 4(c)]. Provided that each atomic layer remains single- 
domain, the net magnetisation vector integrated normal to the plane 
of the junction is independent of position and can be cancelled by the 
appropriate applied field. Therefore, the maximum critical current 
should be independent of the net magnetic moment and so cannot 
account for the dependence reported here. 
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As noted in the introduction, it is also important to consider the 
possibihty that the relative alignment of the Ni layers controls a 
singlet supercurrent. For example, in Ref. 41 it was shown that for 
a Nb-Py-Cu-Co-Nb junction in which an applied field could be 
used to control the relative orientations of the magnetisations of 
the magnetic layers, a much larger critical current was obtained 
when the Py (Ni8oFe2o) and Co layers were antiparallel. A similar 
result was obtained in a related experiment with Fe-Cr-Fe barriers 
in which the Cr spacer layer thickness controlled the relative mag- 
netisation directions of the Fe layers^°; theoretical fits demon- 
strated that this result was consistent with the conventional singlet 
proximity effect. 

In the experiments reported here, the relative alignment of the Ni 
layers is altered by the magnetic history. However, it is obvious from 
Fig. 5 that the configuration in which the Ni layers are antiparallel 
(0 = 180°) has a large net moment (mainly originating from the Gd) 
which corresponds to the minimum rather than the maximum crit- 
ical current. In fact the configuration which corresponds to the max- 
imum critical current for dcd > 1.1 nm is closer to parallel alignment 
of the Ni layers. 

In Fig. 3(e) we have plotted the dependence of the average normal- 
ised supercurrent on the Gd thickness. There are several important 
points to be made about this figure. Firstly, for both the maximum 
and minimum values of the supercurrent, there is a strong reduction 
with increasing Gd thickness, reflecting the short mean free path of 
Gd which translates into a short coherence length for both singlet 
and triplet pairs. Nevertheless, the supercurrent for the 1 .9 nm thick- 
ness data is higher than that for 1.1 nm and comparable with the 
pure Ni sample with no Gd. 

The conclusion is therefore that the significant enhancement in 
the maximum critical current with a reduction in AH that we see in 
our results cannot be due to an incomplete flux cancellation or due to 
a conventional singlet proximity effect. Nevertheless, the experi- 
mental results demonstrate that the junction supercurrent is directly 
dependent on the magnetic structure in the Gd layer. Our analysis 
shows that the supercurrent increases with the angle between the Ni 
layers and hence depends on the helicity of the in-plane DW in the 
Gd layer. Theoretically, in such S-F-F'-F-S junctions, spin-aligned 
triplet pairs form via a spin-mixing effect^'^^. This means that the 
supercurrent cannot be understood purely in terms of simple singlet 
Cooper pair phase-shifts. We argue that the magnetic structure in Gd 
is a Bloch DW, and Eschrig et al in Ref 33 have shown theoretically 
that in the diffusive limit there is an inter- conversion between singlet 
and triplet pairs at the F/DW interfaces (i.e., equivalent to the Ni/Gd 
interfaces in our experiment). The generation of triplet pairs neces- 
sarily means that the Josephson effect is enhanced and so the junc- 
tion critical current increases. Our results are consistent with this 
picture and the increase in the junction supercurrent with increasing 
DW-angle in Gd implies an increasing efficiency in the conversion 
between singlet and triplet pairs. 

Methods 

Growth of multilayer thin-films. Nb-Ni-Gd-Ni-Nb thin-films were grown using 
d.c.-magnetron sputtering in an Ar plasma at 1.5 Pa onto 10 X 4 mm^ and 4 X 
4 mm^ unheated single crystal silicon substrates with a 250-nm-thick oxide layer on 
the surface. Prior to film growth and during the deposition, the walls of the main 
chamber of the system were cooled via a liquid nitrogen jacket to lower the system's 
base pressure to below 10"*^ Pa (verified using an in-situ residual gas analyser). 
Substrates rested on a circular table that rotated below stationary sputtering targets. 
All targets were pre- sputtered until a constant and stable voltage was achieved 
(typically between 10 and 15 minutes). The film growth rate was controlled by the 
target power and the speed in which the substrates passed below the stationary targets. 
Deposition rates and film thicknesses were pre-calibrated using atomic force 
microscopy to measure pre-deposited step-edges. 

Nanopillar device fabrication. Films were patterned by standard optical lithography 
and Ar-ion milling to produce a series of 4-|j,m-wide and 30-|j,m-long wires, each 
connected to four larger area contact pads to allow for four-point measurements to be 
performed. Within each wire a nanopillar Josephson junction was processed using a 
focused Ga-ion beam microscope technique. We have described this process in detail 



elsewhere in Refs. 43 and 44. The junction areas were approximately square and in the 
300 X 300 nm^ to 1200 X 1200 nm^ range. 
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